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ABSTRACT 


In modern yacht design geometric modeling is regarded 
to be directly related to the hydrodynamic performance 
of the shape of the hull and its appending elements - 
usually the keel, often with winglets, and the rudder. 
While the traditional way of shape design - i.e., draw- 
ing, model building, tank testing, modifying ... - is both 
time consuming and expensive, a complementing ap- 
proach shall be discussed within this paper. The ap- 
proach is called hydrodynamic modeling since it tightly 
combines the hydrodynamic analysis and the geomet- 
ric modeling in the design process. It is based on ad- 
vanced Computational Fluid Dynamics (CFD) methods 
for flow field analysis and unique parametric modeling 
techniques for shape generation. 

The geometry of a yacht is entirely described via im- 
portant form parameters as discussed in detail by the au- 
thors at the 1999 CSYS. The canoe body of the yacht is 
modeled from a small set of longitudinal curves which 
provide all parameters needed for sectional design. The 
longitudinal curves themselves being created via form 
parameters, a fully parametric description of the hull is 
achieved which allows to create and modify the geom- 
etry in a highly sophisticated manner. The fairness of 
the shapes is an intrinsic part of the form generation 
procedure. Apart from the canoe body the keel repre- 
sents the most pronounced hydrodynamic design ele- 
ment, dominating lift and righting moment of a yacht 
but also causing a non-negligible resistance component 
called induced drag. Keel, bulb and winglets are also 
specified in terms of form parameters. 

An application of hydrodynamic modeling is given 
for an IACC-yacht. Formal optimization can be suc- 
cessfully employed to identify improved and, eventu- 
ally, optimal configurations. A reasonably small set of 
parameters (free variables) was selected and systemati- 
cally varied making use of a fully-automatic optimiza- 
tion scheme. Two optimization examples are presented 
in order to demonstrate the potential of the approach: (a) 


the optimization of a keel-bulb- winglet configuration so 
as to find a minimum drag solution for a given sideforce 
and (b) the optimization of the bare hull with respect to 
wave resistance. 

The examples can be regarded as representative for 
both racing and touring yachts with draft restrictions and 
illustrate the methodology of hydrodynamic modeling. 

NOMENCLATURE 

B max maximum beam at deck level 

Cb block coefficient 

Cp prismatic coefficient 

£2 fairness criteria 

F„ Froude number 

IACC International America’s Cup Class 
Lpp length between perpendiculars 

VPP Velocity Prediction Program 

x Bma, longitude of maximum beam 

xcb longitudinal center of buoyancy 

Tmax maximum draft 


INTRODUCTION 

Designing a yacht, in particular its hull geometry and 
appendages, is a process of creativity, skill, experience 
and art - independent of whether the naval architect 
chooses to express his or her ideas by means of a tra- 
ditionally drawn lines plan or whether the designer de- 
cides to apply a computer aided design (CAD) system to 
create a product model. Full benefit can be gained from 
the latter when an integrated process of modeling and 
analysis is established in which design variations can be 
evoked and assessed efficiently. 

In geometric modeling and particularly in yacht de- 
sign many CAD systems are now built on an outstand- 
ing mathematical curve and surface representation tech- 
nique known as B-splines. Originating in free-form de- 
sign, the underlying methodology of most of these sys- 
tems is the interactive shape generation where points - 



e.g. the vertices of the B-spline’s defining polygon or 
polyhedron - are positioned in three-dimensional space. 
Achieving the desired form generally is a laborious un- 
dertaking since the results need to be suitably fair while 
specific constraints have to be taken into account, e.g. 
the displacement or the length of a water line important 
to the rating rules under consideration. Then the pro- 
cess of manual vertex manipulation becomes rather te- 
dious and systematic modifications in order to improve 
the shapes with regard to their hydrodynamics become 
inapt. 

Instead of interactively handling the lowest entities 
of the underlying mathematical model (i.e., the vertices) 
a different approach has been pursued which is aimed at 
expressing the geometry in terms of high level descrip- 
tors for the intended shapes (i.e., form parameters). 

Following the stage of shape creation, the design 
may be analyzed for its various characteristics. The hy- 
drodynamics being of supreme importance to racing 
yachts, a state-of-the-art system of Computational Fluid 
Dynamics can be used to examine the performance. 
Modem flow codes have reached the maturity to rank 
different designs in respect to resistance and lift (i.e., 
side force). A potential flow code with good response 
time was therefore applied to numerically analyze the 
flow about the hull and a keel-bulb-winglet configura- 
tion. 

The potential of linking the two stages of geometric 
modeling and hydrodynamic analysis tightly together, 
becomes apparent when utilizing an integrated envi- 
ronment in which modeling, analysis, evaluation and 
modification can be repeated systematically within short 
turn-around time. 

Within this paper the parametric approach will be 
discussed and examples will be shown for an IACC 
yacht, see section on geometric modeling. The hydrody- 
namic optimization of an IACC canoe body and its keel- 
bulb-winglet configuration by means of a formal strat- 
egy and a fully-automatic process will be presented, see 
section on optimization. A full optimization in the light 
of a velocity prediction program (VPP) is discussed. 


GEOMETRIC MODELING 

In parametric modeling design ideas are usually ex- 
pressed by descriptors that imply higher level informa- 
tion about the object to be created. Often, relationships 
and possible dependencies between entities are consid- 
ered. In addition, the descriptors may also represent 
complex features that the product is to assume. When 
modeling geometry, the descriptors are called form pa- 
rameters, three types of which can be distinguished: 

Differential form parameters like tangents and curva- 


tures (e.g. angle of entrance of the design water- 
line), 

Positional form parameters like points to be interpo- 
lated (e.g. breadth of the waterline at the transom), 

Integral form parameters like area, volume and cen- 
troid information (e.g. center of flotation). 

Well-defined parameters facilitate the modeling process 
since the designer can focus his or her attention on the 
outcome rather than on the input, assuming that the 
form generation procedure automatically brings about 
the specified geometry by itself. 

In the subsections to come, first the parametric mod- 
eling of (bare) hulls shall be briefly reviewed as intro- 
duced by Harries and Abt (1999b) and Harries 
(1998). New features will be presented so as to cope 
with additional constraints originating from class rules. 
Following this, the parametric modeling of appendages 
shall be outlined as needed for hydrodynamic optimiza- 
tion. 

Hull 

In the novel parametric approach to the design of sailing 
yachts by Harries and Abt (1999b) the process of 
modeling surfaces of complex geometry is based on lay- 
ing out a set of cross-sectional curves and, subsequently, 
generating a surface by means of lofting (Letcher, 
1981) or skinning (Woodward, 1986, 1988). Follow- 
ing the classic naval architect’s technique of describing 
a ship’s geometry in terms of longitudinal curves (i.e., 
basic curves) from which design sections are derived, 
the design sections (i.e., the cross-sectional curves) de- 
fine the interpolating surfaces and determine the shape 
of the envisioned hull. 

A new system called FRIENDSHIP Modeler} has 
been introduced which is based completely on paramet- 
ric design principles. One of the key features of the 

1 Form parameter oRIENteD SHIP Modeler http://www. 
frieftdsHip- systems . com 
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Figure 1 : Levels of the hydrodynamic modeling sys- 
tem 




Layer 

Objective function 

Constraints 

Design/Free variable 

2 

Hydrodynamic 

optimization 

Performance 

Feasible domains of the 
design variables 

Properties defined by the 
constraints of layer 1 e.g. 
Cp, xcb 

1 

Parametric hull 
generation with implicit 
measurement constraints 
and physical properties 
from form parameter 
defined basic curves 

Fairness of basic curves 
and sectional fairness 

Displacement, xcb , 
interpolation of the deck, 
design waterline , 
centerplane, tangents, 
tangent plane, shape 
modification curves, 
measurement marks 

Direct geometrical 
properties of the basic 
curves e.g. B max , T max , 
Lpp, XB max , xp max , B stern , 
Bbow 

IT 

Parametric B-spline 
generation 

£2 Fairness 

Interpolation, enclosed 
area, centroid position, 
tangential properties, 
curvature 

Vertex coordinates, 
vector sizes 


Table 1 : Direct parametric modeling 


Layer 

Objective function 

Constraints 

Design/Free variable 

3 

Hydrodynamic 

optimization 

Performance 

Feasible domains of the 
design variables 

Properties defined by the 
constraints of layer 2 
(e.g. Cp, xcb ) and free 
variables of layer 1 

~Y~ 

Geometric optimization 

Global and local fairness 
of the hull 

Displacement, formula 
constraints, form 
parameters e.g. Cp, Cb, 
xcb , lateral area, 
waterplane area, center 
of flotation, convexity 

A subset of the free 
variables from layer 1 

— 

Parametric hull 
generation with implicit 
measurement constraints 
and physical properties 
from form parameter 
defined basic curves 

Fairness of basic curves 
and sectional fairness 

Displacement, xcb , 
Interpolation of the deck, 
design waterline , 
centerplane, tangents, 
tangent plane, shape 
modification curves, 
measurement mark 

Direct geometrical 
properties of the basic 
curves e.g. B max , T rnax , 
Lpp, xp max , xp max , B stern , 
Bbow 

IT 

Parametric B-spline 
generation 

£2 Fairness 

Interpolation, enclosed 
area, centroid position, 
tangential properties, 
curvature 

Vertex coordinates, 
vector sizes 


Table 2: Advanced parametric modeling 


approach is the generation of B-spline curves and sur- 
faces by means of variational calculus. Instead of in- 
teractively manipulating the B-spline’s control points, 
the (free) vertices are computed from a geometric op- 
timization which employs fairness criteria as measures 
of merit and captures the specified form parameters as 
equality constraints. 

Modeling a hull thus becomes the task of selecting 
the form parameters to be taken into account and assign- 
ing suitable values to them. This can be done by evalu- 
ating an existing design and remodeling it or, alterna- 
tively, specifying a set of form parameters from scratch. 
As soon as an initial shape is produced changes can be 


systematically brought about by varying one or several 
parameters. 

Form parameters can be individually addressed and 
changed. Nevertheless, the interdependency of form pa- 
rameters needs to be considered. For instance, pushing 
the center of flotation aft while pulling the center of 
buoyancy forward can only be accommodated within 
subtle limits unless non-yacht like shapes are intended. 

Within hydrodynamic optimization, the direct use 
of physical properties like displacement and center of 
buoyancy can be successfully employed, see Harries 
and Abt (1999a). Variations can be evoked efficiently 
but the initial guess should be reasonably close to where 





Figure 2: Parametrically designed IACC yacht with circular sections 



Figure 3 : Parametrically designed IACC yacht with trapezoidal sections 









the optimal shape is expected. Sometimes, however, it 
is beneficial that an initial design is changed consider- 
ably and more than initially assumed. (This might be 
due to the lack of experience when developing an en- 
tirely new project.) Consequently, giving the modeling 
process more freedom will greatly assist in finding an 
optimum - which might then even lie outside the naval 
architect’s conventional experience. 

Furthermore, the resulting shapes may sometimes 
not suit the requirements imposed by class rules. Typ- 
ically, when optimizing for wave resistance on a down- 
wind course at medium Froude numbers bulbous-like 
sections in the vicinity of the bow appear, contradicting 
convexity constraints. (This is due to a favorable redis- 
tribution of displacement - at least for heavy displace- 
ment yachts where dynamic lift does not play an impor- 
tant role.) Considering class rules at the early stage of 
shape generation already establishes a tangible advan- 
tage. 

A special branch of the FRIENDSHIP Modeler - 
called IACC-Friend 2 - has been derived by the authors 
extending the direct modeling approach to comply with 
the IAC-Class rules (IACC, 1997). 

In order to meet the hull’s convexity requirement an 
additional optimization layer has been introduced in the 
design process. Tab. 1 shows the three layers which are 
used in the direct design mode where the sectional area 
curve forms an integral part of the input to the modeling, 
see also Harries and Abt (1999b). Tab. 2 presents 
an advanced approach where large shape variations can 
be realized while class requirements are simultaneously 
fulfilled. An additional layer is introduced to balance the 
shape changes due to a hydrodynamic optimization with 
the restrictions given by the rule. 

For instance, the longitudinal position of the max- 
imum breadth of the design waterline is a typical and 
effective design variable for a hydrodynamic optimiza- 
tion. The longitude of the maximum beam of the deck 
should then be utilized as a free variable at a level where 
formula constraints are accommodated - i.e., at layer 2 
in Tab. 2 - such that the hull maintains its convexity. At 
this level typical integral form parameters which are se- 
lected as design variables are implemented as equality 
constraints. Positional and differential form parameters 
are generally treated at layer 1 but may be passed to ei- 
ther level 2 or 3 depending on the design problem at 
hand. At the lowest level - i.e., at layer 0 - the B-splines 
are computed according to the input received from lev- 
els 1 to 3. 

Because of the underlying fairness criteria employed 
to determine the B-spline curves, the shapes created 
from the parametric approach are intrinsically fair. The 
fairness criteria (see Harries and Abt, 1999b) are 

-International Americas Cup Class Form parameter oRIENteD 


minimized in an optimization, favoring circular sections 
- i.e., natural shapes of roundish character. For bet- 
ter shape control an additional parameter has been de- 
vised which allows to conveniently design sections with 
straight or straightened parts typical of trapezoidal sec- 
tions. This new form parameter is the vertical position 
which a transverse B-spline curve has to interpolate at a 
predefined curve parameter. 

The definition of the form parameters is shown in 
Tab. 3. 


Appendages 

Similar to the direct design mode for modeling the canoe 
body of a yacht, the keel fin and bulb can be parametri- 
cally described and modeled. While the fin and winglet 
may be described readily from excellent wing section 
data - for instance via scaling, blending and merging in 
a longitudinal sweep operation - the bulb generally is 
a free-form object. A bulb’s volume and its distribution 
for example should be accurately defined by means of 
a sectional area curve, enabling the designer to specify 
the mass and center of gravity. 

An advanced feature of this approach is the method 
of modifying natural shapes, i.e., shapes that originate 
from the B-spline optimization when neglecting cen- 
troid information. The centroid of the section is mod- 
ified relative to an unconstrained optimization. This 
means: After computing a section the highest and lowest 
centroid position is determined from a vertex transfor- 
mation in which a zero curvature condition is applied 
at the upper and lower ends of each B-spline curve, 
respectively. From this the extreme centroid locations 
are calculated and mapped onto an unified parameter 
space. Subsequently, each design section is computed 
anew with the desired centroid location retrieved from 
the basic curve which defines the centroid modification 
along the bulb axis. The basic curves defining the con- 
tour, sectional area curve and the centroid modifier are 


Figure 6: Basic curves of an IACC keel bulb 




Figure 7: IACC keel bulb with natural vertical cen- 
troid (neutral centroid modifier) 



Figure 8: IACC keel bulb with modified centroid - low 
to neutral 



Figure 9: IACC keel bulb with modified centroid - 
high via low to neutral 

displayed in Fig. 6 for an example bulb. Naturally, those 
basic curves are also determined via form parameters. 
The bulbs depicted in Fig. 7 to Fig. 9 exactly feature 
the same weight and longitudinal center of gravity but 
originate from changes in the centroid modifier. 

For a hydrodynamic optimization the form parame- 
ters of all four basic curves can be readily applied. The 
volume and center of gravity of the bulb are usually to 
be kept constant while the tangents of the sectional area 



Figure 10: Optimization Process 

curve may be varied. Also, the contours are subject to 
possible change. Of course, from the set of bulb form 
parameters any suitable subset can be selected. 


OPTIMIZATION 

In the proposed parametric design method a yacht’s ge- 
ometry is created in terms of its direct properties as ex- 
pressed by its form parameters. The hull is determined 
from optimizing the fairness criteria and, consequently, 
the generated shapes 

• accurately meet all desired properties and 

• intrinsically acquire excellent fairness. 

Both features are key prerequisites for the optimiza- 
tion of a yacht’s most important indirect properties, i.e., 
its various hydrodynamic qualities like resistance, sea- 
keeping and lift-drag ratio. 

Figure 10 displays the synthesis model for the for- 
mal hydrodynamic optimization of yachts. The synthe- 
sis model comprises four stages: 

Model of form generation: Parametric design via the 
IACC-Friend Modeler. 

Model of hydrodynamic analysis: CFD simulation 
by means of a state-of-the-art flow code. 







Model of design evaluation: Integral lift and drag 
forces. 

Model of optimization: Non-linear programming 
techniques. 

Optimizations take place at two levels. In an outer loop 
important form parameters are systematically varied so 
as to improve a suitable hydrodynamics criterion, e.g. 
the drag for a given sideforce at a predefined boat speed. 
In an inner loop, the geometry is optimized as discussed 
in the previous subsections. 

The optimization scheme in the outer loop is based 
on a conjugate gradient method as described by Press 
ET AL. (1988). The algorithm compromises two steps 
which are alternately repeated until convergence. In the 
first step, the gradient of the measure of merit is com- 
puted with respect to the free variables at a base point. 
In the second step, a promising search direction is iden- 
tified and a one-dimensional optimization is undertaken, 
setting out from the base point into the direction of im- 
provement. Here, the Golden Section search method is 
employed. 

The optimum point found along the search line is 
then used as a new base point and the procedure starts 
anew: the current gradient is computed and a new search 
direction is determined for the next line search. Instead 
of simply using the gradient at the current base point - 
as it would be done in the method of steepest descent - 
the new direction is computed from the current and pre- 
ceeding gradients, promising improved performance in 
long and narrow valleys of the search space. 

In order to demonstrate the feasibility of this ap- 
proach two cases have been studied: 

• A simplified keel-bulb-winglet configuration 
where the position as well as the sweep and 
dihedral angles of the winglets have been varied in 
order to find the configuration with the minimum 
induced drag for a given sideforce. 

• The canoe body where various form parameters 
have been modified so as to identify the minimum 
wave resistance of the yacht sailing upright at a 
moderate Froude number. 

The potential flow module of the CFD system 
SHIPFLOW by Larsson (1997) was used for the numer- 
ical simulations. A reliability study was undertaken for 
sailing yachts previous to the CFD based optimizations, 
see Piller (2000). 

Induced Resistance 

A topological model was created to set up a keel-bulb- 
winglet configuration as shown in Fig. 1 1 which pro- 
vides a finite set of design variables. For simplicity the 



Figure 1 1 : Keel-bulb-winglet configuration used for 
analysis of induced drag (neglecting keel flap) 


fin and winglets were generated from NACA 631-012 
sections (Abbott and von Doenhoff, 1958) while 
the bulb was modeled as described above. Focusing on 
induced drag a double body model was considered to be 
sufficient so as to reduce the computational effort. 

A similar configuration was presented by Larsson 
( 1999), demonstrating that the SHIPFLOW system can be 
successfully applied to the hydrodynamic design of a 
keel-bulb-winglet configuration. 

Keel, bulb and winglets were discretized using 400, 
2000 and 260 panels, respectively, the source strength 
for each panel being defined by higher order distribu- 
tion. In order to fulfil the Kutta-Jukowski condition ad- 
ditional strip groups were added at the trailing edges of 
keel and winglets. Extra strips were introduced to move 
the tip-vortex to the center of the bulb. For details on 
the method see Janson (1997). Fig. 13 displays a typi- 
cal result of these calculations showing velocity vectors 
and pressure contours. 

Calculations have been carried out for three differ- 
ent leeway angles (0.0°, 2.5 ° and 5.0°) and a constant 
heel angle of 22.0 °. The value of induced drag and lift 
was derived from Treffetz plane analysis. The induced 
drag associated with an exemplary sideforce of 25 kN 
was computed via a non-linear interpolation of the re- 
sults for the different leeway angles. 

Only the longitudinal position of the winglets and 
their dihedral and sweep angle were chosen as free vari- 
ables within the optimization run (see Fig. 11) to keep 
the number of computations reasonable. Fig. 12 shows 
the history of the optimization process: within 10 evalu- 
ations the induced resistance could be reduced by ap- 
proximately 2 % - the main change steming from the 
longitudinal position of the winglets. 



Keel Optimization History 



Figure 12: Optimization history for the induced resistance of a keel-bulb- winglet configuration 



Figure 13: Perspective view of a simplified keel-bulb- 
winglet configuration featuring velocity vectors and 
pressure contours from a potential flow calculation 


Wave Resistance 

A second set of calculations was carried out for the bare 
hull of an IACC yacht sailing with neither heel nor lee- 
way. Tab. 3 shows the input file for the parametric mod- 


eler lACC-Friend which completely describes the gen- 
erated hull as depicted in Fig. 5. 

The wave resistance is calculated by SHIPFLOW, in- 
cluding the non-linear boundary condition at free sur- 
face (Fig. 14). For an initial solution approximately 
10 iterations were needed. The modified configurations 
were restarted from a previous solution and converged 
within a few iterations. 

The wave resistance being calculated from pressure 



Figure 14: Wave contours from a non-linear potential 
flow calculation for an IACC yacht 


Hull Optimization History for Fn=0.307 



Figure 1 5 : Optimization history for the canoe body without heel and leeway at F n = 0.307 


integration was regarded as a sensible measure of merit. 
The calculation were carried out for a Froude number of 
F n = 0.307. 

The optimization history is plotted in Fig. 15 for a 
variation of four free variables: 

• longitudinal position of maximum breadth of the 
design waterline 

• bow beam 

• longitudinal position of maximum draft 

• angle of the deck line at the stern 

The wave resistance of the initial configuration is com- 
puted to be 564 N. A tangible reduction of more than 
25 % is then accomplished within the first few calcula- 
tions. Within the following iterations improvements in 
wave resistance can be found as well, however, within 
the range of only a few percent. 

Further variations have been undertaken for different 
starting values of the free variables. These optimizations 
resulted in similar trends. 

Naturally, it should be kept in mind that the achiev- 
able improvement always depends on the quality of the 
initial system. The better the original design the less 
potential for changes. Improving a good design is thus 
much more challenging than gaining on a less sophisti- 
cated initial shape. 


Finally, it needs be pointed out that the accom- 
plished improvements can only be as good as the com- 
putational model employed to analyze the flow field. 
One should therefore not expect the full advantage pre- 
dicted and experimental validation studies are required 
to prove the validity of the hydrodynamic optimization 
based on CFD calculations. 


Yacht Optimization - Outlook 

The performance of a yacht’s hull and its appendages 
can only be fully appreciated when undertaking a com- 
plete velocity prediction (e.g. Hochkirch, 2000), bal- 
ancing the hydrodynamic and aerodynamic forces and 
moments for the entire spectrum of courses and wind ve- 
locities, see Fig. 16. Furthermore, the optimization of a 
yacht under race conditions implies many more compu- 
tations and scenarios than shown here. For instance sea- 
keeping and manoeuvering need to be considered, too. 
Nevertheless, it has not been an attempt of this paper to 
tackle this formidable task but to highlight techniques 
that greatly help in the design process and, eventually, 
facilitate a more comprehensive optimization. 

For an IACC yacht the design evaluation becomes 
a challenging task in itself since a (probabilistic) mea- 
sure of merit ought to be considered. For given wind and 
sea conditions the race time has to be integrated from 
a VPP yielding the amount of seconds needed for the 



Figure 16: Complete optimization process for racing 
yachts 

course. The likelyhood of a specific match race scenario 
in comparison to others will have to be taken into ac- 
count. However, an extraordinary amount of computa- 
tions will be required to determine this ultimate measure 
of merit and to optimize for it. 

CONCLUSIONS 

The design philosophy of hydrodynamic modeling has 
been presented. The approach is based on closely relat- 
ing 

• sophisticated geometric modeling techniques and 

• advanced numerical flow field analysis. 

A synthesis model for hydrodynamic optimization was 
applied which comprises the four stages of form gen- 
eration, fluid dynamic analysis, design evaluation and 
formal optimization. 

Unique parametric modeling techniques have been 
discussed. Particular emphasis has been given to the 


fair design of the canoe body of sailing yachts. An ad- 
vanced parametric description of both the hull and its 
appendages is regarded as the key to successful for- 
mal optimization. Fully-automatic optimization allows 
to study a variety of shapes and, eventually, to improve 
the design with respect to selected measures of merit. 

The parametric design of an IACC yacht as well as 
exemplary optimizations of the yacht’s bare hull and 
keel-bulb-winglet configuration have been shown. The 
parametrically designed shapes feature excellent fair- 
ness. Shape control is accomplished by means flexible 
sets of parameters (the high level descriptors of the de- 
sign ideas). Promising improvements could be achieved, 
demonstrating the potential of hydrodynamic modeling 
of sailing yachts. 

It must be emphasized that the hydrodynamic mod- 
eling approach presented cannot replace the experience 
and skill of the yacht designer. The naval architect still 
is the key figure to decide which parameters to choose 
and vary and to assess the validity of the final design. 
Nevertheless, the approach enables to concentrate on the 
design task instead of the modeling problem. 
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